Ab initio many-body calculations of nucleon scattering on He, 7 Li, 'Be, C and b O 
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We combine a recently developed ab initio many-body approach capable of describing simultane- 
ously both bound and scattering states, the ab initio NCSM/RGM, with an importance truncation 
scheme for the cluster eigenstate basis and demostrate its applicability to nuclei with mass numbers 
as high as 17. Using soft similarity renormalization group evolved chiral nucleon-nucleon interac- 
tions, we first calculate nucleon- He phase shifts, cross sections and analyzing power. Next, we 
investigate nucleon scattering on 7 Li, 7 Be, 12 C and ls O in coupled-channel NCSM/RGM calcula- 
tions that include low-lying excited states of these nuclei. We check the convergence of phase shifts 
with the basis size and study A = 8, 13, and 17 bound and unbound states. Our calculations predict 
low-lying resonances in 8 Li and 8 B that have not been experimentally clearly identified yet. We are 
able to reproduce reasonably well the structure of the A = 13 low lying states. However, we find 
that A = 17 states cannot be described without an improved treatment of 16 one-particle-one-hole 
excitations and a clustering. 



I. INTRODUCTION 

Nuclei are quantum many-body systems with both 
bound and unbound states. A realistic ab initio de- 
scription of light nuclei with predictive power must have 
the capability to describe both classes of states within 
a unified framework. Over the past decade, significant 
progress has been made in our understanding of the 
properties of the bound states of light nuclei starting 
from realistic nucleon-nucleon (NN) interactions, see e.g. 
Ref. [l| and references therein, and more recently also 
from NN plus three- nucleon (NNN) interactions @tJ|- 
The solution of the nuclear many-body problem becomes 
more complex when scattering or nuclear reactions are 
considered. For A = 3 and 4 nucleon systems, the Fad- 
deev and Faddeev-Yakubovsky Q as well as the hyper- 
spherical harmonics (HH) [7| or the Alt, Grassberger and 
Sandhas (AGS) [8| methods are applicable and success- 
ful. However, ab initio calculations for unbound states 
and scattering processes involving more than four nu- 
cleons in total are quite challenging. The first ab initio 
many-body neutron- 4 He scattering calculations were per- 
formed within the Green's Function Monte Carlo method 
using the Argonne NN potential and the Illinois NNN 
interaction [9(. Also, resonances in He isotopes were in- 
vestigated within the coupled-cluster method using the 
Gamow basis (lOj . 

In a new development, we have recently combined 
the ab initio no-core shell model (NCSM) [ll| and the 
resonating-group method (RGM) }12 - ll7l |. into a new 
many-body approach pH E3 {ab initio NCSM/RGM) 
capable of treating bound and scattering states of light 
nuclei in a unified formalism, starting from fundamen- 
tal inter-nucleon interactions. The NCSM is an ab ini- 
tio approach to the microscopic calculation of ground 
and low-lying excited states of light nuclei with realis- 
tic two- and, in general, three-nucleon forces. The RGM 
is a microscopic cluster technique based on the use of A- 
nucleon Hamiltonians, with fully anti-symmetric many- 



body wave functions built assuming that the nucleons are 
grouped into clusters. Although most of its applications 
are based on the use of binary-cluster wave functions, 
the RGM can be formulated for three (and, in principle, 
even more) clusters in relative motion |13| . The use of the 
harmonic oscillator (HO) basis in the NCSM results in 
an incorrect description of the wave-function asymptotic 
and a lack of coupling to the continuum. By combining 
the NCSM with the RGM, we complement the ability 
of the RGM to deal with scattering and reactions with 
the use of realistic interactions, and a consistent ab ini- 
tio description of the nucleon clusters, achieved via the 
NCSM. Presently the NCSM/RGM approach has been 
formulated for processes involving binary-cluster systems 
only. However, extensions of the approach to include 
three-body cluster channels are feasible, also in view 
of recent developments on the treatment of both three- 
body bound and continuum states (see, e.g ., Refs. [20 - 
|24J|). As described in detail in Refs. [Hi DjI, the ab 
initio NCSM/RGM approach has been already applied 



to study the n- H, 



me, 



10 



Be, and p- 3,4 He scat- 



tering processes, and address the parity inversion of the 
11 Be ground state, using realistic NN potentials. In that 
work, we demonstrated convergence of the approach with 
increasing basis size in the case of the A = 4 and A = 5 
scattering. The n- 10 Be calculations were, on the other 
hand, perfomed only in a limited basis due to the com- 
putational complexity of the NCSM calculations of the 
10 Be eigenstates. 

It is the purpose of the present paper to expand the ap- 
plicability of the NCSM/RGM beyond the lightest nuclei 
by using sufficiently large NMl HO excitations to guar- 
antee convergence of the calculation with the HO basis 
expansion of both the cluster wave functions and the lo- 
calized RGM integration kernels. The use of large NMl 
values is now feasible due to the recent introduction of 
the importance truncated (IT) NCSM scheme [25l[2r|. 
It turns out that many of the basis states used in the 
NCSM calculations are irrelevant for the description of 



any particular eigenstate, e.g., the ground state or a set 
of low-lying states. Therefore, if one were able to identify 
the important basis states beforehand, one could reduce 
the dimension of the matrix eigenvalue problem with- 
out losing predictive power. This can be done using an 
importance truncation scheme based on many-body per- 
turbation theory |25j |. 

We make use of the IT NCSM wave functions for 
the cluster eigenstates, in particular the eigenstates of 
the target nucleus of the binary nucleon-nucleus system, 
and calculate the one- and two-body densities that are 
then used to obtain the NCSM/RGM integration ker- 
nels. We benchmark the IT approach in basis sizes ac- 
cessible by the full calculation and apply it within still 
larger basis sizes until convergence is reached for tar- 
get nuclei as heavy as 12 C or 16 0. In this study, we 
employ a similarity renormalization group (SRG) [271 - 
H evolved chiral N 3 LO NN potential J3c| (SRG-N 3 LO) 
that is soft enough to allow us reach convergence within 
about 14 — 16M1 HO excitations in the basis expansion. 

In Sect. HI we briefly overview the NCSM/RGM for- 
malism and present for the first time the IT NCSM 
scheme that includes both ground and low-lying excited 
states in the set of reference states. Next, we present 
scattering calculation results for the n- 4 He and p- 4 He 
systems in Sect. IIII1 In particular, we compare the cal- 
culated phase shifts to an R-matrix analysis of exper- 
imental data and, further, calculated differential cross 
sections and analyzing powers in the energy range 6- 
19 MeV to the corresponding experimental data. Neu- 
tron elastic and inelastic scattering on 7 Li and proton 
elastic and inelastic scattering on 7 Be are investigated 
in Sect. IIVI We present phase shifts, cross sections and 
scattering lengths. We predict resonances in 8 Li and 8 Be 
that have not been clearly identified in experiments yet. 
In Sect.[V] we discuss nucleon- 12 C results for both bound 
and unbound states of 13 C and 13 N, obtained includ- 
ing two 12 C bound states, the ground and the first 2 + 
state, in the NCSM/RGM coupled-channel calculations. 
In Sect. I VII we present results for the nucleon- 16 sys- 
tem. In the NCSM/RGM coupled-channel calculations, 
we take into account the 16 ground state and up to the 
lowest three 16 negative-parity states. Conclusions are 
given in Sect. I VIII 
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In the above expression, [A—aail^Ti) and laazT^'lz) 
are the internal (antisymmetric) wave functions of the 
first and second cluster, containing A— a and a nucleons 
(a<A), respectively. They are characterized by angular 
momentum quantum numbers I\ and I<x coupled together 
to form channel spin s. For their parity, isospin and ad- 
ditional quantum numbers we use, respectively, the no- 
tations 7Ti,Tj, and ai, with i = 1,2. The cluster centers 
of mass are separated by the relative coordinate 
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where {fi,i — 1,2, ••• , A] are the A single-particle co- 
ordinates. The channel states ((T|) have relative angu- 
lar momentum £. It is convenient to group all rele- 
vant quantum numbers into a cumulative index v — 
\A.—aa\I* x T\\ a«2/ 2 ' r2 T2; s£}. 

The former basis states can be used to expand the 
many-body wave function according to 
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As the basis states ([l) are not anti-symmetric under ex- 
change of nucleons belonging to different clusters, in or- 
der to preserve the Pauli principle one has to introduce 
the appropriate inter-cluster anti-symmetrizer, schemat- 
ically A v = 



- . /{A-ay.a 



— r — X)p(~) P -P' where the sum runs 
over all possible permutations P that can be carried out 
among nucleons pertaining to different clusters, and p 
is the number of interchanges characterizing them. The 
coefficients of the expansion ([3]) are the relative-motion 
wave functions g'l T (r), which represent the only un- 
knowns of the problem. To determine them one has 
to solve the non-local integro-differential coupled-channel 
equations 



II. FORMALISM 



A. NCSM/RGM 
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(4) 
where the two integration kernels, the Hamiltonian ker- 
nel, 



The ab inito NCSM/RGM approach was introduced in 
Ref . [l8| with details of the formalism given in Ref . [19J ■ 
Here we give a brief overview of the main points. 

In the present paper, we limit ourselves to a two-cluster 
RGM, which is based on binary-cluster channel states of 
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contain all the nuclear structure and anti-symmetrization 
properties of the problem. In particular, the non-locality 
of the kernels is a direct consequence of the exchanges of 
nucleons between the clusters. We have used the notation 
E and H to denote the total energy in the center-of-mass 
frame, and the intrinsic A-nucleon microscopic Hamilto- 
nian, respectively. 

The formalism presented above is combined with the 
ab initio NCSM in two steps: 

First, we note that the Hamiltonian can be written as 



1 2 ■ ■ ■ A-2 A A-l 



H 



T re i(r) + V rcl + V c (r) + H {A _ a) + H {a) , (7) 



where H( A _ a j and H/ a -\ are the (A— a)- and a-nucleon 
intrinsic Hamiltonians, respectively, T re i(r) is the rela- 
tive kinetic energy and V rc i is the sum of all interac- 
tions between nucleons belonging to different clusters af- 
ter subtraction of the average Coulomb interaction be- 
tween them, explicitly singled out in the term Vc{r) = 
Z\ v Z2 V e 2 jr (Z\ v and Z 2v being the charge numbers of 
the clusters in channel v). We use identical realistic 
potentials in both the cluster's Hamiltonians and inter- 
cluster interaction V re i ■ Accordingly, \A— aail^Ti) and 
laa^i^ 2 ^} are obtained by diagonalizing Hi A _ a \ and 
-ff( a ), respectively, in the model space spanned by the 
NCSM N max htt HO basis. Note that in the present pa- 
per we use soft SRC evolved NN potentials. Therefore, 
there is no need to derive any further effective interac- 
tion tailored to the model space truncation as with these 
soft interactions our calculations converge in the model 
spaces we are able to reach. 

Second, we replace the delta functions in the localized 
parts of the Hamiltonian ([5]) and the norm (j6|) kernels 
with their representation in the HO model space. We 
use identical HO frequency as for the cluster eigenstate 
wave functions and a consistent model space size (iV max ). 
We emphasize that this replacement is performed only 
for the localized parts of the kernels. The diagonal parts 
coming from the identity operator in the antisymmetriz- 
ers, the kinetic term and the average Coulomb potential 
are treated exactly. 

In this paper, we apply the NCSM/RGM formalism 
in the single-nucleon projectile basis, i.e., for binary- 
cluster channel states ([lj with a — 1 (with channel index 
v = {A— lail^Tx] l||; s£}). As an illustration, let's 
discuss in more detail the norm kernel that is rather sim- 
ple in this basis: 
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We can easily recognize a direct term, in which initial and 
final state are identical (corresponding to diagram (a) of 
Fig.[T]), and a many-body correction due to the exchange 
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(c) 



(d) 



(e) 



FIG. 1: Diagrammatic representation of: (a) "direct" and (b) 
"exchange" components of the norm kernel; (c and d) "direct" 
and (e) "exchange" components of the potential kernel. The 
group of circled black lines represents the target cluster, a 
state of A— 1 nucleons. The separate red line represents the 
projectile, a single nucleon. Bottom and upper part of the 
diagram represent initial and final states, respectively. 



part of the inter-cluster anti-symmetrizer (corresponding 
to diagram (6) of Fig. [T]). We note that in calculating 
the matrix elements of the exchange operator Pa—i,a we 
replaced the delta function of Eq. (fT|) with its represen- 
tation in the HO model space as discussed above. This is 
appropriate as the transposition Pa-i,a operator acting 
on the traget wave function is short-to-medium range. 
On the contrary, the (5-function coming from the identity 
is treated exactly. The presence of the inter-cluster anti- 
symmetrizer affects also the Hamiltonian kernel, and, in 
particular, the matrix elements of the interaction. For a 
TV TV potential one obtains a direct term involving inter- 
action and exchange of two nucleons only (see diagrams 
(c) and (d) of Fig. [IJ , and an exchange term involving 
three- nucleons. Diagram (e) of Fig. Q] describes this lat- 
ter term, in which the last nucleon is exchanged with one 
of the nucleons of the first cluster, and interacts with 
yet another nucleon. For more details on the integration 
kernels in the single-nucleon projectile basis we refer the 
readers to Ref. fl9j j . 

Being translationally-invariant quantities, the norm 
and Hamiltonian kernels can be "naturally" derived 
working within the NCSM Jacobi-coordinate basis. How- 
ever, by introducing Slater-determinant channel states of 
the type 
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in which the eigenstates of the (A— a)-nucleon fragment 
are obtained in the SD basis (while the second cluster 
is still a NCSM Jacobi-coordinate eigenstate), it can be 
easily demonstrated that translationally invariant matrix 
elements can be extracted from those calculated in the 



SD basis of Eq. (fit)]) by inverting the following expression: 
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Here Ot.i. represents any scalar and parity-conserving 
translational-invariant operator (Ot.i. = A, AHA, etc.). 
We exploited both Jacobi-coordinate and SD channel 
states to verify our results. The use of the SD basis is 
computationally advantageous and allows us to explore 
reactions involving p-shell nuclei, as done in the present 
work. In order to calculate the parts of the integration 
kernels depicted in Fig.Q](b), (c) and (d), all information 
that we need from the SD basis calculation are one-body 
densities of the target eigenstates. For the (e) part of the 
integration kernel in Fig. [TJ we need two-body densities 
of the target eigenstates obtained in the SD basis. 

Due to the presence of the norm kernel .A/' I ^ J / r (r / ,r), 
Eq. (QJ does not represent a system of multichan- 
nel Schrodinger equations, and g J v T (r) do not repre- 
sent Schrodinger wave functions. The short-range non- 
orthogonality, induced by the non-identical permutations 
in the inter-cluster anti-symmetrizers, can be removed by 
introducing normalized Schrodinger wave functions 



xf T {r) 



E J dyy 2 Nhr,y) 



gf r (y) 
y 



(12) 



where A/" 5 is the square root of the norm kernel, and ap- 
plying the inverse-square root of the norm kernel, A/" - 3 ( 
to both left and right-hand side of the square brackets 
in E q. ffl . This procedure, explained in more detail in 
Ref. [19| , leads to a system of multichannel Schrodinger 
equations 



[T re i(r) + V c (r) - (E - e£ Tx - E$ T2 )]^M 
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I 7Ti Ti 

where E a \ ' is the energy eigenvalue of the i-ih cluster 
(i = 1,2), and W^,J'{r' ,r) is the overall non-local po- 
tential between the two clusters, which depends on the 
channel of relative motion, while it does not depend on 
the energy. These are the equations that we finally solve 
to obtain both our scattering and bound-state results. 



B. Importance truncated NCSM with excited 
states 

The primary limitation for the range of applicability 
of the NCSM in terms of particle number A and model 
spaces size iV max results from the factorial growth of the 
dimension of the iV max M7 space. Except for light iso- 
topes, it is hardly possible to obtain a converged result 
using a 'bare' Hamiltonian within the N max HQ spaces 
that are computationally tractable. 

At this point the importance truncation offers a solu- 
tion. The importance truncation in connection with the 
NCSM was introduced in Ref. [25[ and discussed in de- 
tail in Ref. 26] . In the following we summarize a few key 
features of the IT-NCSM and generalize the approach to 
the simultaneous description of excited states. 

The motivation for the importance truncation results 
from the observation that the expansion of any particular 
eigenstate of the Hamiltonian in a full m-scheme NCSM 
space typically contains a large number of basis states 
with extremely small or vanishing amplitudes. The am- 
plitudes define an adaptive truncation criterion, which 
takes into account the properties of the Hamiltonian and 
the structure of the eigenstate under consideration. If 
those amplitudes were known — at least approximately — 
before actually solving the eigenvalue problem, one could 
reduce the model space to the most relevant basis states 
by imposing a threshold on the amplitude. The ampli- 
tude of a particular basis state |$„) in the expansion of a 
specific eigenstate can be estimated using first-order mul- 
ticonfigurational perturbation theory. In order to set up 
a perturbation series we need an initial approximation 
of the target state, the so-called reference state |^ re f). 
In practice this reference state will be a superposition of 
basis states |$„) £ M^i from a reference space M rc f. 



hM= E ^ ref) l $ M> 
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The reference state and the amplitudes C^ v are typi- 
cally extracted from a previous NCSM calculation. Based 
on |\& rc f) as unperturbed state, we can evaluate the first- 
order perturbative correction to the target state resulting 
from basis states |$„) £ M rc f- Their first-order ampli- 
tude defines the so-called importance measure 
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The energy denominator e„ — e re f in a M0ller-Plesset- 
type partitioning is simply given by the unperturbed 
harmonic-oscillator excitation energy of the basis state 
|$„) (see Ref. for details). 

Imposing an importance threshold K m ; n , we construct 
an importance truncated model space including all basis 
states with importance measure \k v \ > K m in- Since the 
importance measure is zero for all basis states that differ 
from all of the states in A4 rc { by more than a two-particle- 
two-hole excitation, we have to embed the construction of 



the importance truncated space into an iterative update 
cycle. After constructing the importance truncated space 
and solving the eigenvalue problem in that space, we ob- 
tain an improved approximation for the target state that 
defines a reference state for the next iteration. In order 
to accelerate the evaluation of the importance measure 
(fT5)l . we typically do not use the complete eigenstate as 
new reference state, but project it onto a reference space 
spanned by the basis states with \C V \ > C mm , where 
C v are the coefficients resulting from the solution of the 
eigenvalue problem. The second threshold C m - ln will be 
chosen sufficiently small so as not to affect the results for 
a given K m ; n threshold. 

Simple iterative update schemes can be devised for any 
type of full model spaces, as discussed in Refs. [261 [31] . 
Specifically for the N max hfl space of the NCSM, how- 
ever, there is an efficient sequential update scheme lead- 
ing to the IT-NCSM(seq) approach. It is based on the 
fact that all states of an (-/V max + 2)Ml space can be gen- 
erated from the basis states of an iV max Ml space using 
two-particle-two-hole excitations at most. Thus a single 
importance update starting from a reference state in an 
N max _Ml space gives access to all relevant states in an 
(Amax + 2)Ml space. Making use of this property, in the 
IT-NCSM(seq) we start with a full NCSM calculation in, 
e.g., 2MI and use this eigenstate after applying the C m { n 
threshold as reference state for constructing the impor- 
tance truncated AMI space. After solving the eigenvalue 
problem for this importance truncated AMI space we use 
the resulting eigenstate as reference state to construct 
the QMl space, and so on. Thus only one importance 
update is required for each value of -/V max , which makes 
this scheme very efficient computationally. Moreover, in 
the limit of vanishing thresholds, (« m i n , Cmin) — > 0, this 
scheme recovers the full N max hil space at each step of 
the sequence, i.e., the IT-NCSM(seq) would recover the 
full NCSM result. 

Based on this limiting property, we can obtain a nu- 
merical approximation to the full NCSM result by extrap- 
olating the IT-NCSM(seq) observables obtained for a set 
of different importance thresholds K m i n (and in princi- 
ple also C m ; n ) to K m in — > 0. Through this extrapolation, 
the contribution of discarded basis states, i.e. those with 
importance measures \k v \ below the smallest threshold 
considered, is effectively recovered. Because the control 
parameter K min is tied to the physical structure of the 
eigenstate, we observe a smooth threshold dependence 
for all observables, which allows for a robust threshold 
extrapolation. In the case of the energy we can improve 
the quality of the extrapolation further by considering a 
perturbative second-order estimate for the energy of the 
excluded basis states. While setting up the importance 
truncated space, all second-order energy contributions 
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value. By construction this correction goes to zero in the 
limit K m ; n — > 0. We use this additional information for a 
constrained simultaneous extrapolation of the energy to 
vanishing threshold with and without perturbative cor- 
rection for the excluded states as described in detail in 
Ref. [26|. 

The whole concept can be generalized to the descrip- 
tion of excited states. For the present application in con- 
nection with the NCSM/RGM, we aim at an importance 
truncated model space that is equally well suited for the 
description of the lowest M eigenstates of the Hamilto- 
nian for given parity and total angular momentum pro- 
jection. Instead of using a single reference state, we em- 
ploy different reference states l^lf ), with m = 1, ...,M, 
for each of the M target states. For each reference state 
we define a separate importance measure Ky following 
Eq. (fT5|) . A basis state |$„) is included in the impor- 
tance truncated space if at least one of the importance 
measures |«i, m | is above the threshold K m i n , i.e., if it is 
relevant for the description of at least one of the M tar- 
get states it will be included. Because the different eigen- 
states are typically dominated by different basis states, 
the dimension of the importance truncated space grows 
linearly with M. 

In the IT-NCSM(seq) scheme we start with a full 
NCSM calculation in 2MI and use the lowest M eigen- 



states as initial reference states |\I> 
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for the discarded states with \k v | < K min are summed up 
to provide a correction A exc i(K m ; n ) to the energy eigen- 



corresponding importance measures k v the importance 
truncated AMI space is constructed and the lowest M 
eigenvectors in this space serve as new reference states 
(after application of the C m i n threshold) for the construc- 
tion of the 6M7 space, and so on. From a sequence of IT- 
NCSM(seq) calculations we obtain a set of M eigenvec- 
tors for each value of ./V max which can be used to evaluate 
other observables. 

By default we compute the expectation values of J 2 
and T 2 as well as the expectation values of H lnt and H cm . 
Indeed, since we use an importance truncated space in 
the m-scheme without explicit angular momentum pro- 
jection, the eigenstates are not guaranteed to have good 
angular momentum and isospin. We therefore monitor 
the expectation values of J 2 and T 2 and find values which 
typically differ by less then 10 -3 from the exact quan- 
tum numbers. As in the full NCSM we separate spurious 
center-of-mass (CM) excitations from the physical spec- 
trum by adding a Lawson term (3H crn to the translation- 
ally invariant intrinsic Hamiltonian Hi nt (with the typi- 
cal choice j3 = 10). The use of this modified Hamiltonian 
provides at the same time a diagnostic for potential CM 
contaminations of the intrinsic states induced by the im- 
portance truncation. As discussed in Refs. [26|, [32|], the 
independence of the intrinsic energies (-Hint) on /3 and the 
smallness of (H cm ) demonstate that the IT-NCSM(seq) 
solutions are free of CM contaminations. 

Eventually, the wave functions obtained in the IT- 
NCSM(seq) together with the threshold extrapolated in- 
trinsic energies form the input for the NCSM/RGM cal- 



dilations discussed in the following. 



III. NUCLEON- 4 HE SCATTERING 

The purpose of the nuclcon- 4 He calculations presented 
in this paper is two-fold. First, we want to check the 
predictive power of the SRG evolved chiral interaction in 
the A = 5 system, where a lot of experimental scatter- 
ing data exist and where our calculations can be easily 
converged with respect to the size of the basis expansion. 
Second, we want to benchmark the importance trunca- 
tion scheme with the full-space calculations all the way 
up to very large iV max ft$7 spaces. 

The first ab initio A = 5 scattering calculations was 
reported in Ref. Q. The n-a low-lying J 71 ' = 3/2 _ and 
l/2~ P-wave resonances as well as the 1/2+ S'-wave non- 
resonant scattering below 5 MeV cm. energy were ob- 
tained using the AV18 NN potential with and without 
the three-nucleon force, chosen to be either the Urbana 
IX or the Illinois-2 model. The results of these Green's 
function Monte Carlo (GFMC) calculations revealed sen- 
sitivity to the inter-nucleon interaction, and in particular 
to the strength of the spin-orbit force. 

Soon after, the development of the ab initio 
NCSM/RGM approach allowed us to calculate both n- 
and (for the first time) p-a scattering phase shifts for en- 
ergies up to the inelastic threshold [lg, [l9j , using several 
realistic N N p otentials, including the chiral N 3 LO [301 ) . 
the Vi owfc [H| and the CD-Bonn [H NN potentials. 
Nucleon-a scattering provides one of the best-case sce- 
narios for the application of the NCSM/RGM approach. 
This process is characterized by a single open channel 
up to the <i+ 3 H threshold, which is fairly high in energy. 
In addition, the low-lying resonances of the 4 He nucleus 
are narrow enough to be reasonably reproduced diago- 
nalizing the four-body Hamiltonian in the NCSM model 
space. In the present work we include the first excited 
state of 4 He, the + state, as a closed channel in our 
NCSM/RGM basis space. 
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FIG. 2: (Color online) Dependence of the n- 4 He phase shifts 
on the size of the HO basis expansion of the He wave func- 
tions and the localized parts of the integration kernels. The 
4 He ground state and the first + excited states were in- 
cluded. The SRG-N 3 LO NN potential with a cutoff of 2.02 
fm - and the HO frequency Ml = 20 MeV were used. 



that obtained using the V\ ow k NN potentatial in our er- 
lier study (compare the present Fig. [2] to the left panel 
of Fig. 13 in Ref. [H]). 



Benchmark of Importance- Truncated 
calculations 



A. Convergence with the size of the HO basis 
expansion 

We performed extensive nucleon- 4 He calculations with 
the SRG-N 3 LO NN potential with a cutoff of 2.02 fm- 1 
to check convergence of our NCSM/RGM calculations. 
In Fig. [2J we present n- 4 He phase shift results for the S- 
and P- waves obtained using an HO basis expansion up to 
A max = 17 for for the localized parts of the NCSM/RGM 
integration kernels and for the 4 He ground- and the first- 
excited + wave functions (since these states have posi- 
tive parity, the N m&x — 1 expansion is in fact used for the 
4 He eigenstates). As seen in the figure, the phase-shift 
convergence is excellent. In particular, the N max = 17 
and the N max — 15 curves lie on top of each other. The 
convergence rate demonstrated here is quite similar to 



As shown in the previous subsection, for the A = 5 
system we are able to reach complete convergence with 
4 He wave functions obtained within full, non-truncated, 
NCSM calculations. We can, therefore, test the perfor- 
mace of the IT-NCSM scheme in this system all the way 
up to very large iV max values and see how well the IT- 
NCSM scheme reproduces the completetly converged re- 
sults. It should be noted that for heavier A = 8, 13 and 
A = 17 systems investigated later, full, non-truncated 
NCSM calculations for the A = 7 (A = 12, 16) target 
nuclei are feasible only up to -/V max = 10 (JV max = 8). It 
is, therefore, desirable and important to benchmark the 
IT-NCSM calculations in a lighter system like A = 5 in 
iV mM > 10 calculations. 

In Fig. [31 we compare n- 4 He phase shifts calculated 
within the NCSM/RGM with 4 He wave functions ob- 
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FIG. 3: (Color online) Calculated n- 4 He 5- and P-wave phase 
shifts. Results obtained with 4 He wave functions from full 
NCSM (solid lines) and IT-NCSM (dashed lines) calculations 
are compared. The SRG-N LO NN potential with a cutoff of 
2.02 fm _1 , the iV max = 17 basis space and the HO frequency 
Ml = 20 MeV were used. See text for details on the IT-NCSM 
calculation. 
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FIG. 4: (Color online) Calculated n— 4 He (left panels) and 
p— 4 He (right panels) compared to the R-matrix analysis of 
experimental data [35|]. The NCSM/RGM calculations that 
included the 4 He ground state and the + excited state were 
done using the SRG-N 3 LO NN potential with a cutoff of 2.02 
fin" 1 . The HO frequency Ml = 20 MeV and iV max = 17 basis 
space were employed. 



C. Comparison with experimental data 



tained in a full N max = 16 NCSM calculation and those 
obtained using He wave functions obtained within an 
iV max = 16 IT-NCSM calculation. The agreement of the 
two sets of phase shifts is excellent. It should be noted 
that the dimension of the full N max = 16 4 He NCSM ba- 
sis is 6344119. The dimension of the IT-NCSM basis used 
here to calculate the 4 He wave functions was just 992578, 
more than a factor of six smaller. Truncation parameters 
K min = 10 -5 and C min = 2 x 10~ 4 were used. The ground 
state energy from the full NCSM calculation is —28.224 
MeV. The n min — > extrapolated ground state energy 
from the IT-NCSM calculation is -28.217(5) MeV with 
a difference from the full result less than 10 keV. The 
excited + energy obtained in the full NCSM calcu- 
lation was 21.58 MeV. The corresponding extrapolated 
IT-NCSM result was 21.4(1) MeV. The slightly lower ac- 
curacy of the excited state reproduction in the IT-NCSM 
calculation is manifested in a very small deviation of the 
S-wave phase shift at energies above 12 MeV (less than 
1 degree at 16 MeV). It should be noted that the excited 
0+0 state is not bound. Consequently, it is challenging 
to reproduce the excited state as well as the ground state 
in a importance-truncated calculation. It should be also 
pointed out that unlike for the energies, no phase shift 
extrapolation was performed. The needed one- and two- 
body densities were calculated from the wave functions 
obtained in the IT-NCSM calculation with the truncation 
parameters described above. The excellent agreement of 
the full and the IT-NCSM phase shifts demonstrates that 
no extrapolation was actually necessary. Obviously, we 
can check the dependence of observables like phase shifts 
on the K m j n and C m i n and perform an extrapolation to 
vanishing values of these parameters if needed. 



Our calculated n- 4 He and p- 4 He phase shifts are com- 
pared to those obtained from an i?-matrix analysis of 
N— 4 He experimental data [35|] in Fig. [H The agree- 
ment is quite reasonable for the S'-wave, £>-wave and 
2 Pi/2-wave. The 2 Ps/2 resonance is positioned at higher 
energy in the calculation and the corresponding phase 
shifts are underestimated with respect to the P-matrix 
results, although the disagreement becomes less and less 
pronounced starting at about 8 MeV. While the inclusion 
of negative-parity excited states of the a— particle would 
likely increase somewhat the 2 Pz/2 phase shifts [18l[l9l]. 
the observed difference is largely due to a reduction in 
spin-orbit strength caused by the neglect of the three- 
nucleon interaction in our calculations. The importance 
of the three-nucleon force in reproducing the P-matrix 
2 P 3 /2 phase shifts was demonstrated in the GFMC n- 
4 He calculations of Ref. [9j. Overall, the present results 
obtained with the SRG-N 3 LO NN interaction agree bet- 
ter with experiment than our earlier calculations [let I19J 
with the Mowfc, N 3 LO and CD-Bonn NN potentials. The 
only exception is the S'-wave phase shift which is best 
described using the CD-Bonn NN potential. The larger 
spin-orbit strength of the employed SRG-N 3 LO potential 
with respect to N 3 LO itself is the likely responsible for 
the improved agreement. 

As our calculated phase shifts agree with the exper- 
imental ones reasonably well above the center-of-mass 
energy of 8 MeV, we expect a similar behavior for cross 
section and analysing power in that energy range. This 
is indeed the case as shown in Fig. where the calcu- 
lated differential cross section and analyzing power are 
compared to experimental data from Karlsruhe [36j with 
polarized neutrons of E n —17 MeV laboratory energy. 
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FIG. 5: (Color online) Calculated n— 4 He differential cross section for neutron laboratory energy of E n — 17 MeV (a), and 
analyzing power for E n — 17 (&), 15 (c) and 19 MeV (d) compared to experimental data from Ref. 36]. The NCSM/RGM 
results include the 4 He ground state and first + excited state and were obtained using the SRG-N 3 LO AW potential with a 
cutoff of 2.02 fm" 1 , for an HO frequency Ml = 20 MeV and basis space size iV max = 17. 
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FIG. 6: (Color online) Calculated p— 4 He differential cross section (top panels) and analyzing power (bottom panels) for proton 
laboratory energies of E p - 5.95, 7.89,9.89 and 11.99 MeV compared to experimental data from Ref. [37|. The NCSM/RGM 
results include the 4 He ground state and first + excited state and were obtained using the SRG-N 3 LO AW potential with a 
cutoff of 2.02 fm~ , for an HO frequency hfl — 20 MeV and basis space size iV max = 17. 
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FIG. 7: (Color online) Calculated p— 4 He differential cross section (top panels) and analyzing power (bottom panels) for proton 
laboratory energies E p — 12, 14.32 and 17 MeV compared to experimental data from Schwandt et al. [37(, Brokman [3g|, Dodder 
et al. [39| . and Hardekopf et al. [401 ]. The NCSM/RGM results include the 4 He ground state and first + excited state and 
were obtained using the SRG-N 3 LO AW potential with a cutoff of 2.02 fori -1 , for an HO frequency h£l = 20 MeV and basis 
space size iV ma x = 17. 



For the cross section experimental data see also refer- 
ences in [36j . The cross section is reproduced remarkably 
well at all angles and the analysing power is in reason- 
able agreement with the data, particularly at backward 
angles. The same quality of agreement can be found for 
all energies far from the low- lying resonances, as shown 
in the right panel of Fig. [S] for the analysing power at 
E n = 15 MeV and 19 MeV. 



A better display of the dependence of our calculated 
cross section and analysing power upon the incident nu- 
cleon energy is provided by Fig. [5J where the p— 4 Hc re- 
sults for these observables are compared to the data of 
Ref. [37J at the proton laboratory energies of E p = 5.95, 
7.89, 9.89, and 11.99 MeV. As expected from the be- 
havior of the phase shifts described earlier, for energies 
relatively close to the resonance region we find a rather 
poor agreement with experiment, particularly noticeable 
in the analysing power overall and in the cross section at 
backward angles. However, starting at about 10 MeV, 
the agreement improves substantially and data are once 
again reproduced in a quite satisfactory way at higher 
energies, as shown in Fig. [71 where the NCSM/RGM 
p— 4 He results are compared to various experimental data 



sets [37H40| in the energy range E p ~ 12 - 17 MeV. 



IV. NEUTRON- 7 LI AND PROTON- 7 BE 
SCATTERING 

The 7 Be(p,7) 8 B capture reaction plays a very impor- 
tant role in nuclear astrophysics as it serves as an input 
for understanding the solar neutrino flux [4jJ. While the 
experimental determination of the neutrino flux from 8 B 
has an accuracy of about 9% [42[ , the theoretical predic- 
tions have uncertainties of the order of 20% [43, |44j . The 
theoretical neutrino flux depends on the 7 Be(p,7) 8 B S- 
factor. Significant experimental and theoretical effort has 
been devoted to studying this reaction. The S-factor ex- 
trapolation to astrophysically relevant energies depends 
among other things on the scattering lengths of the pro- 
ton scattering on 7 Be. Experimental determination of 
these lengths was performed recently [45[ with precision 
of the order of 30%. The proton- 7 Be elastic scattering 
was also investigated in Ref. [46] . To benchmark the the- 
oretical calculations used for S-factor extrapolations, an 
investigation of the mirror capture reaction, 7 Li(n,7) 8 Li, 
as well as the n+ 7 Li scattering is important. For exam- 
ple, the n+ 7 Li scattering lengths are known with a higher 
accuracy |47| . 

The first applications of the NCSM approach to the de- 
scription of the 7 Be(p,7) 8 B capture reaction [48| required 
a phenomenological correction of the asympotic behav- 
ior of the overlap functions and, further, the scattering 
p+ 7 Be wave function was calculated from a phenomeno- 
logical potential model. The present investigation within 
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FIG. 8: (Color online) 7 Li ground-state and the l/2~ and 
7/2~ excited state energy dependence on the model-space 
size iV max , obtained within the importance-truncated NCSM 
(solid lines), using the SRG-N 3 LO NN potential with a cutoff 
of 2.02 fin" 1 . The HO frequency Wl = 20 MeV was employed. 
The full-space NCSM results are shown by dashed lines. 



the ab initio NCSM/RGM approach paves the way for a 
complete first principles calculation of this capture reac- 
tion. Here, we limit ourselves to scattering calculations 
and postpone the capture reaction calculations to a forth- 
coming paper. 

Our current limit on the unrestricted NCSM calcula- 
tions for 7 Li and 7 Be is N max = 10. To improve the con- 
vergence of our scattering calculations, we utilize wave 
functions obtained within the IT-NCSM. In that scheme, 
we are able to reach V max = 18 model spaces and calcu- 
late both ground as well as low-lying excited states. This 
is demonstrated in FigM With the SRG-N 3 LO NN po- 
tential with A = 2.02 fm _1 employed in the present study 
we reach convergence already around iV max = 12 — 14. 
Also, as seen in the figure, the aggreement between the 
unrestricted NCSM and the IT-NCSM is perfect up to 
the highest accessible unrestricted space, V max = 10. 



A. n- 7 Li 

The NCSM/RGM coupled-channel calculations per- 
formed for the A — 8 system include the 7 Li ( 7 Be) ground 
state, the first excited l/2~ state as well as the second 
excited 7/2~ state. It is essential to include the 7/2~ 
state in order to reproduce the low-lying 3 + resonance 
in 8 Li and 8 B. Using these three states, we are able to 
reach model spaces up to -V max = 12, which is sufficient 
concerning the HO basis expansion convergence as can 
be judged from Fig. [5] The coupled channel calculation 
described above gives two bound states for the n- 7 Li sys- 
tem, a 2 + corresponding to the experimentally observed 
8 Li ground state, bound by 2.03 MeV [49], and a 1 + cor- 
responding to the 8 Li first excited state at E x = 0.98 
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FIG. 9: (Color online) P-wave diagonal phase shifts of the 
n- 7 Li elastic scattering (top panel), elastic 7 Li(n,n) 7 Li cross 
section (middle panel), and inelastic 7 Li(n,n') 7 Li(l/2~) cross 
section (bottom panel). The NCSM/RGM calculation that 
included the 7 Li ground state and the l/2~ and 7/2~ excited 
states were done using the SRG-N 3 LO NN potential with a 
cutoff of 2.02 fin -1 . Wave functions from IT-NCSM calcula- 
tions in the iV max = 12 basis and the HO frequency of Ml = 20 
MeV were employed. Experimental data are from Ref. [521 ]. 



MeV, bound by 1.05 MeV [49]. The calculated states 
are bound by 1.16 MeV and 0.17 MeV, respectively, i.e. 
less than in experiment. This is in part due to the 
fact that higher excited states of 7 Li were omitted. In 
Fig. [31 we present our results for the diagonal P-wave 
phase shifts of the n+ 7 Li elastic scattering as well as 
the elastic 7 Li(n,n) 7 Li and inelastic 7 Li(n,n') 7 Li(l/2~) 
cross sections. At low energies, we can identify four res- 
onances two of which can be associated with the experi- 
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mentally known 8 Li states: 3 + at E x = 2.255 MeV and 
1+ at E x = 3.21 MeV J49|. The other two resonances 0+ 
and 2 + are not present in the 8 Li evaluation of Ref. [49j ■ 
They do appear in many theoretical calculations includ- 
ing the GFMC p|, NCSM [H] and recoil-corrected con- 
tinuum shell model (RCCSM) [j^. The 0+ resonance 
also appears in the GCM calculations of Ref. 51]. Con- 
tributions of different resonances to the cross sections 
can be deduced from Fig. [9j The elastic cross section is 
dominated by the 3 + resonance with some contributions 
from the 2 + rcsonace at higher energy. The inelastic 
cross section shows a peak just above the threshold due 
to the + resonance and also a contribution from the 
1 + resonance. The appearance of a + peak just above 
threshold of the 7 Li(n,n') 7 Li(l/2~) reaction was also dis- 
cussed in Ref. [501 ( see ^ig- 10 in that paper). The data 
of Ref. [52j seem to rule out a + state so close to the 
threshold. It is known, however, that the position of the 
+ state is sensitive to the strength of the spin-orbit in- 
teraction 0, [48], [50J ■ The three- nucleon interaction, that 
would increase the strenght of the spin-orbit force, was 
not included in our present calculations. Consequently, 
our predicted + state energy is likely underestimated. 
We note that no fit to the experimental treshold was 
done in the present NCSM/RGM calculations. Still, as 
seen in the bottom panel of Fig. [9] the calculated inelas- 
tic cross section is very close to the experimental data 
just above the threshold. 



B. p- 7 Be 

In the mirror system, p- Be, we do not find a bound 
state in the same type of coupled-channel NCSM/RGM 
calculation as described above for n- 7 Li. As seen in 
the top and the middle parts of Fig. [TOl the lowest 2 + 
resonance corresponding to the 8 B ground state lies at 
about 200 keV above the threshold. In experiment, B 
is bound by 137 keV J49|. Our calculated lowest 1 + res- 
onance appears at about 1 MeV. It corresponds to the 
experimental 8 B 1+ state at E x = 0.77 MeV (0.63 MeV 
above the p- 7 Be threshold). This resonance dominates 
the inelastic cross section as seen in the bottom part of 
Fig. [10] The higher lying resonances follow similar pat- 
terns as those found in n- 7 Li (Fig. [9J. Again, we find 
+ and 2 + resonances not included in the recent 8 B eval- 
uation 49]. We note that experimental efforts are now 
under way to find these resonances [40, [53| . We further 
note that our calculated lj states in 8 Li and 8 B appear 
at a significantly higher energies than the corresponding 
lj" states obtained within the microscopic cluster model 
in Ref. 0. 

The middle panel of Fig. [10] demonstrates once again 
the good accuracy of the importance truncated calcula- 
tions for a high NhQ, N max = 10, model space. The IT 
calculation reduced the 7 Be basis from 43.6 million to 
11.9 million in the present case. 

The elastic p- 7 Be scattering was measured at 148° and 
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FIG. 10: (Color online) P-wave diagonal phase shifts of the 
p+ 7 Be elastic scattering (top and middle panel) and inelas- 
tic 7 Be(p,p') 7 Be(l/2~) cross section (bottom panel). The 
NCSM/RGM calculation that included the 7 Be ground state 
and the 1/2" and 7/2" excited states were done using the 
SRG-N 3 LO NN potential with a cutoff of 2.02 fin -1 . Wave 
functions from IT-NCSM calculations in the iV max = 12 basis 
and the HO frequency of ZlO = 20 MeV were employed. In 
the middle panel, the full-space NCSM (solid lines) and the 
IT-NCSM (dashed lines) results in the JV max = 10 basis are 
compared. 



analyzed by the R-matrix approach |46|. Cross section 
calculations within the RCCSM at that angle were then 
published in Ref. [55J and also in Ref. [50j. Further, elas- 
tic and inelastic cross sections at this angle were analyzed 
within the timc-dcpcndcnt approach to the continuum 
shell model (TDCSM) [H. 0ur elastic and inelastic dif- 
ferential cross section results at 148° are presented in 
Fig. HU In the elastic cross section, the first 1 + state is 
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FIG. 11: (Color online) Elastic 7 Be(p,p) 7 Be (top panel) and 
inelastic 7 Be(p,p') 7 Be(l/2~) (bottom panel) differential cross 
section at 6 c . m . = 148° calculated within the NCSM/RGM 
with SRG-N 3 LO NN potential with A = 2.02 fm _1 . 



visible and beyond the minimum of the cross section, we 
can see the dominant peak due to the 3 + state. At higher 
energies, the 2 + state contributes as well. The inelastic 
cross section at 148° has a similar shape as the reaction 
cross section shown in Fig. [TO] The first 1 + state peak 
dominates at low energy with contributions from the + 
and the second 1 + at higher energies. Our findings are 
in line with the RCCSM results. However, we remind 
the reader that there is no fitting in our calculations, all 
results being predictions based on the SRG-N 3 LO NN 
potential. Because of this, the positions of our calcu- 
lated resonances, e.g., 1 + , 3 + do not exactly reproduce 
experiment. We do not include the experimental data 
in the figure as they would be shifted compared to the 
calculated peaks. There are at least two reasons why our 
predictions do not match the experimental resonances 
accurately. First, our nuclear Hamiltonian is incomplete, 
e.g. no three-nucleon interaction is included. Second, we 
omitted higher resonances of 7 Li and 7 Be due to numer- 
ical reasons. Most likely, the omitted resonances would 
produce some shifts in the calculated peaks. Both these 
points can and will be addressed in the future. Still, our 
current results contain the bulk of the physics behind the 
investigated scattering processes. 
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FIG. 12: (Color online) S-wave phase shifts of the n+ 7 Li 
(solid lines) and the p+ 7 Be (dashed lines) elastic scattering. 
The calculations as described in Figs. [9l and 1101 



C. S-wave scattering lengths of n- 7 Li and p- 7 Be 

In Fig. [T^l we present our calculated n- 7 Li and the p- 
7 Be iS-wave phase shifts. We do not find any evidence for 
a 2~ resonance advocated in Ref. [46| and discussed in 
Ref. [53]. The corresponding scattering lengths together 
with the experimental values are given in Table HI With 
the exception of the p- 7 Be aoi, which has a large experi- 
mental uncertainty, our calculated scattering lengths do 
agree with experimental data as to their signs, there are 
however differences in the absolute values. Again, as di- 
cussed above, the results presented here serve only as a 
first step towards the ab initio investigation of the rt- 7 Li 
and the p- 7 Be reactions. Prospects for a realistic calcu- 
lation of the 7 Be(p,7) 8 B capture are excellent. Here we 
found the 8 B unbound by only 200 keV. It is quite pos- 
sible that 8 B will become bound (with the NN potential 
employed here: SRG-N 3 LO with A = 2.02 fin" 1 ) by in- 
cluding more excited states of 7 Be in the coupled-channel 
NCSM/RGM calculations. The effect of higher excited 
states of 7 Be can be, in fact, most efficiently included by 
coupling the presently used NCSM/RGM basis with the 
8 B NCSM eigenstates as outlined in Ref. [H]. Even if 
the 8 B would not be bound or, most likely, the theshold 
energy will not agree with the experiment, we have the 
possibilty to explore a variation of the SRG NN poten- 



7 Li 



Calc. 



Expt. 



7 Be 
Calc. Expt. 



a 01 [fm] 
a 02 [fm] 



+1.23 
-0.61 



+0.87(7) 
-3.63(5) 



-1.2 
-10.2 



25(9) 
-7(3) 



TABLE I: The n- 7 Li and the p- 7 Be S-wave scattering lengths. 
Theoretical values correspond to calculations as described in 
Figs. [?J1 and 1101 Experimental values are from Refs. [la. |47j. 
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tial evolution parameter A and tune this parameter to fit 
the experimental threshold. We note that for any A the 
SRG-evolved NN potential will describe all two-nucleon 
properties as accurately as the original starting NN po- 
tential, here the chiral N 3 LO potential of Ref. [30|. It 
should be noted that by adding the three-nucleon inter- 
action, omitted in the present calculations due to com- 
putational reasons, the need for a fine-tuning should be 
significantly reduced, i.e. the results should become A 
independent. 



V. NUCLEON- 12 C SCATTERING 

For nucleon scattering calculations on 12 C or heav- 
ier targets within the NCSM/RGM, the use of the im- 
portance truncation becomes essential. For 12 C, the 
full-space NCSM calculations are currently limited to 
N max — 8 (although successful runs were already per- 
formed for AT max = 10 on the biggest supercomputers 
with the latest version of the code MFD j59|). This is 
insufficient for reaching or approaching convergence of 
the 12 C NCSM calculations as seen from Fig. [13] and 
even more so of the NCSM/RGM scattering calcula- 
tions. The importance-truncated calculations, on the 
other hand, are feasible up to N max = 18, where con- 
vergence is reached for both the ground state as well ex- 
cited states. Our 12 C calculations are performed with the 
SRG-N 3 LO NN potential with the evolution parameter 
A = 2.66 fin -1 , a higher value (i.e. shorter evolution, 
less soft) than that used for the lighter nuclei. The use 
of a small A results in large overbinding of heavier nuclei 
and a significant underestimation of their radii. As seen 
in Fig. 1131 our converged 12 C binding energy is about 
84.5(8) MeV, smaller than the experimental value of 92 
MeV and, further, the agreement of the full-space and 
importance-truncated results is perfect all the way up to 

Nmsnr = 8. 



O— O Importance-Truncated NCSM 
0- ♦ + Full NCSM 
•— • 2 + Importance-Truncated NCSM 
O- O 2 + Full NCSM 




FIG. 13: (Color online) Ground-state and the first excited 
2 + state energy dependence on the model-space size iV max for 
12 C, obtained within the importance-truncated NCSM, using 
the SRG-N 3 LO NN potential with a cutoff of 2.66 fin -1 . The 
HO frequency Ml — 24 MeV was employed. The calculation 
is variational. No NCSM effective interaction was used. The 
full NCSM results were obtained with the code Antoine 16011. 
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Our low-energy p- 12 C phase shift results are shown 
in Fig. Q3] The comparison of the A^ max = 16 and 
Nmax = 14 results demonstrates good convergence with 
respect to the HO basis expansion. The 12 C ground 
state and the first 2 + state were included in the coupled- 
channels NCSM/RGM equations. We note that we also 
performed a phase shift comparison of the full-space and 
the importance-truncated calculations up to A^ max = 6 
and found a similarly perfect agreement as presented in 
Fig. [3] for n- 4 He. In the present p- 12 C calculations, we 
found a single bound state, l/2~ at -2.98 MeV, corre- 
sponding to the 13 N ground state, bound experimentally 
by 1.94 MeV [6l| . The lowest resonance in our calcula- 
tion is 3/2~, barely visible at 0.25 MeV above threshold. 
In experiment, this resonance is at 1.56 MeV. Our cal- 
culated l/2 + resonance appears at about 1.5 MeV above 



FIG. 14: (Color online) Thep- 12 C eigenphase shifts calculated 
within the NCSM/RGM using the SRG-N 3 LO NN potential 
with a cutoff of 2.66 fni -1 and the HO frequency ft£l = 24 
MeV. Full lines (dotted lines) corespond to results obtained 
in the JV max = 16 (iV max = 14) model space. The ground 
state and the first excited 2 + state of 12 C was included. The 
1 C wave functions were obtained within the IT NCSM. 



treshold (in experiment at 0.42 MeV above threshold) 
and the 5/2 + resonance at about 4.9 MeV (in experi- 
ment at 2.61 MeV). 



B. 



In the mirror system, n- C, our NCSM/RGM calcu- 
lations produce three bound states: l/2~ at -5.34 MeV 
corresponding to the 13 C ground state experimentally 
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FIG. 15: (Color online) The n- 12 C phase shifts calculated FIG - 16: ( Color online ) The analysing power for n- 12 C elas- 
tic scattering below and above the calculated 5/2 resonance. 
Energies are in the center of mass. The calculation as de- 
scribed in Fig. 1151 



within the NCSM/RGM using the SRG-N 3 LO NN poten- 
tial with a cutoff of 2.66 fin" 1 . The HO frequency Ml = 24 
MeV and the model-spaces size of JV max = 16 were used. The 
ground state and the first excited 2 + state of 12 C was in- 
cluded. The 12 C wave functions were obtained within the IT 
NCSM. 



bound by 4.95 MeV with respect to the n- 12 C thresh- 
old, 3/2~, bound by 2.23 MeV (experimentally bound 
by 1.26 MeV), and 1/2+ bound by just 0.03 MeV (ex- 
perimentally bound by 1.86 MeV). In experiment, there 
is also a 5/2 + state bound by 1.09 MeV. Our present 
NCSM/RGM calculations including the lowest 0+ and 
and the lowest 2 + 12 C states do not produce any bound 
5/2+ state. 

Our low-energy n- 12 C diagonal phase shifts are shown 
in Fig. [151 The 5/2+ resonance is found at 2.8 MeV 
(experimenally at 1.92 MeV with respect to the n- 12 C 
threshold). The steep drop of the 1/2+ phase shift is due 
to the presence of the very weakly bound 1/2+ state. 
We note that similarly as in the case of n Be, discussed 
in Ref. 18], we observe a significant decrease of the 1/2+ 
state energy in the n- 12 C NCSM/RGM calculation when 
compared to the standard NCSM calculation for 13 C. We 
were able to make these comparisons in model spaces up 
to N max — 6 where we found this drop to be about 3 
MeV. 

Analysing powers were measured for proton and neu- 



tron scattering on C [62H64J ] and scattering experiments 
on polarized proton target are under way [65j . In Fig. 1161 
we present our calculated analyzing power below and 
above the energy of the 5/2+ resonance. We note that 
our calculated 5/2+ resonance appears at 2.8 MeV in 
the center of mass (experimentally at 1.92 MeV). Be- 
low the resonance, the analyzing power is positive at 
©cm < 90° and negative at 6cm > 90°. At ener- 
gies above the resonance, the analyzing power reverses 
its sign. Similar observations were made in calculations 
performed within the multichannel algebraic scattering 
(MCAS) theory [H, |67|. See in particular Fig. 5 of 
Ref. |3 . 



Our calculated 13 N and 13 C bound-state levels and 
resonances are more spread than the experimental ones. 
This is a consequence of an underestimation of the 12 C ra- 
dius found to be 2.05 fm with the SRG-N 3 LO NN poten- 
tial. To remedy this, one would have to calculate three- 
nucleon interaction terms induced due to the SRG evolu- 
tion. This can be done as described in Ref. [68(. However, 
we still need to further develop the NCSM/RGM formal- 
ism in order to handle three-nucleon interactions in the 
scattering calculations. 



VI. NUCLEON- 16 SCATTERING 

The calculation of nucleon scattering on 16 is the 
most challenging among the systems we investigate in 
this paper. The a clustering plays an important role in 
the structure of 16 0, in particular for the first excited 0+ 
state that is known to be almost impossible to reproduce 
in NCSM or coupled-cluster calculations. Our present 
calculations do not include the a clustering yet. 

As in the case of 12 C, we rely on the importance- 
truncated NCSM calculations for obtaining the 16 wave 
functions as the full- JVmax NCSM calculations are possi- 
ble only up to iV max = 8. In Fig.[T71 we show the ground- 
state convergence within the IT-NCSM and a comparison 
to the full-space results. Again, up to the largest accessi- 
ble model space, the agreement between the importance- 
truncated and the full-space calculations is perfect. 



O 



It is straightforward to converge nucleon- 16 scatter- 
ing calculations within the NCSM/RGM using the HO 
expansion up to A max = 18. Our calculated n- 16 phase 
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FIG. 17: (Color online) Ground-state energy dependence 
on the model-space size -/V max for 16 0, obtained within the 
importance-truncated NCSM, using the SRG-N 3 LO NN po- 
tential with a cutoff of 2.66 fm _1 . The HO frequency M7 = 24 
MeV was employed. The calculation is variational. No NCSM 
effective interaction was used. The full NCSM results were 
obtained with the code Antoine 16011. 
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FIG. 18: (Color online) The n— 16 phase shifts calculated 
within the NCSM/RGM using the SRG-N 3 LO NN potential 
with a cutoff of 2.66 fm^ 1 and the HO frequency Ml = 24 
MeV in the -/V max = 18 model space. The ground state and 
of 16 was included. The le O wave functions were obtained 
within the IT NCSM. 



shifts are shown in Fig. [T5] and the HO-basis expan- 
sion convegence is checked for the S- and the D-w&ve 
in Fig. [TTJ] These calculations included the 16 ground 
state only. We find two bound states, l/2 + at -0.88 MeV 
and 5/2+ at -0.41 MeV with respect to the n- 16 thresh- 
old. In experiment, the 17 ground state is 5/2 + , bound 
by 4.14 MeV, and the l/2 + state is the first excited state 
bound by 3.27 MeV. There are also two additional bound 
states: 1/2 - and 3/2 _ . Those are unbound in our calcu- 
lations. 

Clearly, it is insufficient to consider only the ground 
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FIG. 19: (Color online) Basis size dependence of the n— 16 
phase shifts calculated within the NCSM/RGM using the 
SRG-N 3 LO NN potential with a cutoff of 2.66 fin -1 . The HO 
frequency of Ml = 24 MeV was used. The ,F = l/2+(3/2+) 
channel is shown in the top (bottom) panel. Model space 
sizes up to JV ml i = 18 were considered. The ground state and 
of 16 was included. The 16 wave functions were obtained 
within the IT NCSM. 



state of 16 in the coupled-channel NCSM/RGM scat- 
tering calculations. We, therefore, include in addition 
the three lowest 16 negative parity states: 3 - , l - , and 
2~. Due to computational limitations, in this case we 
used HO basis expansion up to -/V max =13. Comparing 
Fig. (2H1 to Fig. [T51 the lp — lh negative-parity excited 
states of 16 generate negative-parity resonances in 17 0. 
These resonances do appear, however, at much higher 
energy than in experiment. The reason for this is the 
fact that our calculated 16 lp — lh states have too large 
excitation energy. In particular, our calculated 3 _ ex- 
cited state has an excitation energy of 15.99 MeV while 
experimentally it lies at just 6.13 MeV. One reason for 
the discrepancy is the softness of the SRG-N 3 LO NN 
potential we use that results in an overall overbinding 
of the 16 ground state and in an underestimation of 
its radius. Another aspect is the challenging problem of 
the IT-NCSM extrapolations of the independent positive 
and negative-parity state calculations. The uncertain- 
ties of the relative excitation energies are higher than in 
same-parity calculations. On the positive side our calcu- 
lation with the negative-parity states, even though with 
overestimated excitation energies, results in the proper 
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FIG. 20: (Color online) The n- 16 phase shifts calculated 
within the NCSM/RGM using the SRG-N 3 LO AW potential 
with a cutoff of 2.66 fm" 1 and the HO frequency Ml = 24 
MeV in the Ar max = 13 model space. The ground state and 
and the lowest 3~, 1~ and 2~ excited states of 16 were 



included. The 
IT NCSM. 



'O wave functions were obtained within the 



ordering of the 17 bound states. The ground state is 
5/2+ at -1.32 MeV and the l/2 + state gains binding as 
well, appearing at -1.03 MeV. 



O 



We also investigated the p- 16 scattering and 17 F 
states. When the NCSM/RGM calculations are re- 
stricted to the channels involving only the 16 ground 
state, we find a 1/2+ resonance at 1.0 MeV and a 5/2+ 
resonance at 2.2 MeV. These resonances correspond to 
the 17 F 1/2+ first excited state, bound by 0.105 MeV, 
and the 17 F 5/2 + ground state bound by 0.6 MeV with 
respect to the p+ w O threshold. By coupling channels 
involving the lp — lh 16 3~, 1~ and 2~ excited states, 
the calculated l/2 + and 5/2+ states are still unbound 
resonances but their energy moves significantly closer to 
the threshold: the 1/2+ appears at +0.7 MeV and the 
5/2+ at +1.2 MeV. 

The 17 F low-lying states were recently investigated 
within the coupled-cluster approach with the Gamow- 
Hartree-Fock basis [69j[ . In those calculations with the 
N 3 LO NN potential, the 1/2+ state is weakly bound 
while the 5/2+ state remains unbound by about 0.1 MeV. 
Using the SRG evolved interaction, the 5/2+ state be- 
came bound with the decrease of the cutoff A. We note 
that our calculated 16 ground state energy, -139.0(8) 
MeV (Fig. [17]) obtained with the SRG-N 3 LO NN poten- 
tial with A = 2.66 fm™ 1 , compares well with the CCSD 
coupled-cluster 16 calculations: -137.6 MeV with the 
SRG-N 3 LO NN potential with A = 2.8 fm" 1 [70]. The 
differences in the positions of the 1/2+ and the 5/2+ are 
due to deficiencies in our description of the negative par- 



ity lp— lh states, which could be related to the two-body 
Hamiltonian used here as well as the uncertainties of the 
threshold extrapolations for the excitation energies. The 
inclusion of additional 16 excited states would increase 
the absolute energy of our calculated 17 F states. The 
most efficient way to do this is by coupling the presently 
used NCSM/RGM basis with the 17 F NCSM eigenstates 
in as outlined in Ref. [5al. 



VII. CONCLUSIONS 

By combining the importance truncation scheme 
for the cluster eigenstate basis with the ab initio 
NCSM/RGM approach, we were able to perform many- 
body calculations for nucleon scattering on nuclei with 
mass number as high as A = 16. With the soft SRG- 
evolved chiral NN potentials, convergence of the calcula- 
tions with respect to the HO basis expansion of the target 
eigenstates and the localized parts of the NCSM/RGM 
integration kernels can be reached using AT max = 12 — 16. 

We first benchmarked the IT-NCSM results with the 
full-space NCSM results for the A = 5 system. Our 
neutron- 4 He and proton- 4 He calculations compare well 
with an R-matrix analysis of the data in particular at 
energies above 8 MeV, and describe well measured cross 
sections and analysing powers for those energies. 

Our calculations of n- 7 Li and p- 7 Be scattering predict 
low-lying 0+ and 2+ resonances in 8 Li and 8 B that have 
not been experimentally clearly identified yet. We found 
that the prospects of a realistic ab initio calculation of the 
7 Be(p,7) 8 B capture within our approach are very good. 
In the present calculations we found the 8 B unbound by 
only 200 keV. It is quite possible that 8 B will become 
bound (with the NN potential employed here: SRG- 
N 3 LO with A = 2.02 fm -1 ) by including more excited 
states of 7 Be in the coupled-channel NCSM/RGM calcu- 
lations. Even if the 8 B will still not be bound or, most 
likely, the threshold energy will not agree with the exper- 
iment, we have the possibility to explore a variation of 
the SRG AW potential evolution parameter A and tune 
this parameter to fit the experimental threshold. 

The use of the importance-trunacted basis becomes es- 
sential in calculations with 12 C or 16 targets as the full- 
space NCSM calculations are limited to Ar max = 8. Our 
n- 12 C and p- 12 C investigations included 12 C ground and 
the first excited 2+ states. We found a single bound state, 



1/2+ 



13 



N as in experiment. In 3 C, we found three 



bound states with the 5/2+ state still unbound contrary 
to experiment. Our calculated spectrum of A = 13 states 
is more spread than in experiment due to the underesti- 
mation of the 12 C radius, a consequence of the softness 
of the SRG-evolved NN interaction. 

The description of nucleon scattering on 16 within our 
formalism was the most challenging. The a clustering 
that plays an important role in the structure of 16 is 
not yet included in our present calculations. Further, the 
lp— lh 16 excited states are more difficult to treat in the 
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IT-NCSM approach, as the extrapolations of excitation 
energies are done from the independent ground state and 
the negative-parity state calculations. We found a strong 
impact of the lp — \h 16 states on the positions of the 
lowest A — 17 states. For example, correct ordering of 
the 5/2 + and the l/2 + states in 17 was obtained only 
when the lp — \h states were included. 

Overall, we find that the inclusion of additional excited 
states of the target nuclei would be beneficial in all stud- 
ied systems and more significant with the increase of A. 
Coupled-channel NCSM/RGM calculations with many 
excited states of the traget are computationally challeng- 
ing. The most efficient way of including the effects of such 
states is by coupling the presently used NCSM/RGM ba- 
sis, consisting of just a few lowest excited states, with the 
NCSM eigenstates of the composite system as outlined in 
Ref. [5g. Work on this coupling is under way. 

The use of the SRG-evolved NN interaction facilitates 
convergence of the NCSM/RGM calculations with re- 
spect to the HO basis expansion. On the other hand, 
due to the softness of these interactions, radii of heav- 
ier nuclei become underestimated. To remedy this, one 
would have to calculate three-nucleon interaction terms 
induced due to the SRG evolution. This can be done 
as described in Ref. roll. It is essential to further de- 



velop the NCSM/RGM formalism in order to handle 
three-nucleon interactions, both genuine and the SRG- 
cvolution induced, in the scattering calculations. 

In the present paper, we limited ourselves to 
singlc-nucleon projectile scattering. Extensions of the 
NCSM/RGM formalism to include deuteron, 3 H and 3 He 
projectiles are under way. 
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